Dissipative properties of vibrated granular materials 
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We investigate collective dissipative properties of vibrated granular materials by means of molec- 
ular dynamics simulations. Rates of energy losses indicate three different regimes or "phases" in the 
amplitude-frequency plane of the external forcing, namely, solid, convective, and gas-like regimes. 
The behavior of effective damping decrement in the solid regime is glassy. Practical applications are 
dicussed. 
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The dominating approach in the world of vibration 
control and suppression by granular systems has been 
mainly practical Q. Granular motion relaxes rapidly 
once the energy supply is switched off, and dampers can 
efficiently absorb energy released by shocks of external 
forcing. Engineers classify granular dampers as passive 
ones. 

For the case of "granular gases", i.e. particulate sys- 
tems in a state where the mean free path is large as com- 
pared with particle sizes the cooling rate, the dissipation 
rate of the system has been investigated |^ , and applica- 
tions of this work require an analysis of granular gas (hy- 
dro)dynamics in a given experimental setup. Damping 
in dense granular arrangements is a much more difficult 
problem which is mostly studied experimentally. 

In this work, by using molecular dynamics simulations, 
we show that granular systems reveal different damp- 
ing regimes indicating collective dissipation modes. Our 
study of these regimes leads to a "phase diagram" of 
horizontally vibrated granular systems (see Fig. By 
using this diagram along with the presented estimates 
for damping decrements, practitioners may accelerate the 
design and testing procedures. 

In simulations we focus on two dimensional contain- 
ers which are partially filled with granular material and 
shaken horizontally. The motion of the container is si- 
nusoidal, x(t) — Asin{Lot); it mimics practical situations 
where dampers are tested in the vicinity of the eigen- 
modcs of the vibrating mechanism. We study the reac- 
tion of the system to the choice of parameters of shaking 
A and w, keeping all other parameters (size, roughness 
and hardness of particles, filling factor, size and shape of 
the apparatus) fixed 

Our primary objective is the rate of energy dissipation, 
computed by cycle averaging under steady conditions of 
oscillatory motion. Dissipation is obtained by using two 
different ways to ensure consistency of the data: (i) from 
the total power transmitted to the container walls, and 
(ii) from the dissipative work in inter-particle collisions. 
Numerically, both results coincide within a few percent. 

For the molecular dynamics simulations, we use a 



modified soft-particle model by Cundall and Strack [Q: 
Two particles i and j, with radii Ri and Rj and po- 



sition vectors fi and fj, interact if their compression 
= Ri + Rj — \ fi — rj\ is positive. In this case the 
colliding spheres feel forces and Ff^ ^ , in normal 

and shear directions denoted by unit vectors and , 
respectively. 
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and the resulting momenta acting upon the particles are 
Mi = FfjRi, Mj = -FfjRj. The constant a is the char- 
acteristic dissipation rate of the material, Y is the Young 
modulus and v the Poisson ratio. The normal and shear 
friction coefficients, 7„ = aY/{l—h'^) and 75, model dissi- 
pation during particle contact. Eq. (^ takes into account 
that the particles slide upon each other for the case that 
the Coulomb condition Aij^ij^j < \Ffj\ holds, otherwise 



they feel some viscous friction. R^f = R,Rj/ {Ri + Rj) 
is the effective radius, and the effective mass m^f^ is de- 



fined analogously, 
contact is 
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with fli and flj being the angular velocities of the parti- 
cles. 

The values of the coefficients used in simulations are 
y/(l - :/2) ^ 7 5 ^ -^q7^ 7„ = 7 X 102, 75 = 30, At = 0.5. 
Cgs units are implied throughout the paper. 

With the system parameters specified above, depend- 
ing on forcing one may find intensive convection (Fig. |^). 
Convection patterns in horizontally vibrated granular 
material have been recently reported ||^. In our case 
only two rolls could be observed. Different aspect ratios 
or material parameters give different convection patterns, 
with 2, 4 or more convection rolls (not shown here) 
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FIG. 1. Images of the cycle- averaged motion of the system vibrated at / = 20 for different velocities of forcing, 
Alu: a) 1900, b) 1250, c) 600, d) 100, e) 15. 



The velocity profiles in Fig. |^ have been produced by 
averaging particle velocities in the regime of steady oscil- 
lations. We do not intend to discuss the effect of convec- 
tion in horizontally shaken material in detail, although 
we note that the onset of convection in the system is due 
to a critical amplitude of driving velocity {Alo)^^. Fig. || 
shows the maximum absolute value of convective motion 
in the system, i.e. the length of the longest arrow in 
Fig. |l|, over the velocity amplitude Auj for different fre- 
quencies / (in Hz). Below {Auj)^^ « 60 there is slow 
collective motion in the system, but close to this point 
there is a transition into a rapid convective regime (see 
also Fig. 0d). 

To characterize the dissipation in the system we intro- 
duce an effective damping parameter h which is propor- 
tional to the ratio between the averaged dissipated power 
per cycle T and the mean translational kinetic energy of 
the granular system [B[, 



6 = 



T /t Wdiss{t)dt 
2E„ mnjj,vl{t)dt 



(4) 



Here summation is performed over all granular particles 
in the container. In the reference case of linear oscilla- 
tor the damping decrement b equals the inverse of the 



amplitude relaxation time. 
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FIG. 2. Maximum velocity of collective flow V over 
amplitude of shaking velocity Alu, for the different fre- 
quencies analyzed: o : 5; □ : 10; o : 20; > : 40; A : 80 
(Hz). 
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Fig. |3| shows the damping 6 as a function of the effec- 
tive acceleration F = AuP' / g and the amphtude of the 
velocity of vibration, Auj. Different symbols (filled or 
open) display different frequencies. Except for the very 
low frequency range, b scales with the amplitude of the 
velocity of vibration, Au (Fig. ^). A transition from 
the Aw-scaling into another regime takes place at F « 1 
(Fig. ||b). Below this point the effective damping pa- 
rameter becomes less sensitive to F and fluctuates very 
strongly. 
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FIG. 3. Effective damping parameter 6 vs the ampli- 
tude of velocity of vibration, Auj (a); vs the dimensionless 
forcing acceleration parameter, F (b). Different symbols 
refer to different frequencies; < : 2.5; o : 5; □ : 10; o : 20; 
[> : 40; A : 80 (Hz). Open symbols refer to fast cool- 
ing schedule, filled symbols refer to slow cooling (see text 
for details). The solid line in (a) shows a fit using (^) 
and (^). Observe the fit lines with slopes -1 and -1-1 for 
h oc {Auj)~^ and b oc {Aui)'^^m the "liquid" and "gas" 
regimes, resp. 



Fig. ^ reveals three regions, indicative of different dy- 
namic regimes. The local minimum in Fig. ^ is the re- 
gion of spatially organized behavior with well developed 
convection rolls, where the entire granular mass partici- 
pates in collective motion (Fig. |l|c). For higher velocity 
amplitudes Alo^ the dissipation rate increases, as parti- 
cles begin to fly across the box and the ordered structure 
of the rolls starts to disappear (cf. Fig. ||a-b). The mini- 
mum corresponds to a characteristic velocity at which the 
system begins to display the "gas" state. This character- 
istic velocity can be estimated by equating the character- 
istic time of motion in horizontal and vertical directions 
for a particle with velocity of the order of Aw, 

L/Aljj ^ Auo/g . (5) 

For L = 100, this gives Au — 313, which gives an order- 
of-magnitude estimate of the minimum of Fig. |^a. Addi- 
tional analysis reveals that the correct prefactor is close 
to 2 (see Fig. |^ and discussion below). At higher ve- 
locities particles continue to stay airborne, and gravity 
becomes unimportant. Thus, the boundary Auj ~ VLg 
separates the "liquid" and gas-like regimes |p^ . 

To the left from the minimum, with amplitude decrease 
at any frequency, the depth of the rolls diminishes pro- 
gressively until a critical value {Auj)cr ~ 60 is reached. 
At this velocity the rolls vanish completely and we do not 
find organized motion in the system anymore (cf. Figs. ^ 
and|l|d). The system seems to be in a "solid" state. The 
critical value of velocity is manifested in Fig. |^a by a 
change in the damping slope. Given that the accelera- 
tion of shaking here exceeds g, critical velocity in this 
range points to Auj \/Rg, where R is the particle ra- 
dius. This is consistent with our particle sizes 

In the solid regime one has to switch from velocities to 
accelerations. Following Fig. ^ to small accelerations, 
one finds a change in the behavior at F « 1, i.e. when 
the maximal acceleration of shaking becomes comparable 
with gravity g. For F ^ 1 the curves are smooth, but for 
F < 1 suddenly the data become very noisy. At this point 
the grains form a glassy "solid" phase. We note that our 
understanding of the transition at F ^ 1 is different from 
the solid-fluid transition reported recently |l^ as long as 
Auj ^ \/Rg- At higher velocities our results lead to the 
same conclusions as the results reported in Ref. 

To reiterate, the transition point found in Ref. |pl| is 
given by Auj'^ — g in our notation. It is unclear whether 
the non-glassy solid phase at F > 1 can be identified by 
the method of granular temperature used in Ref. [[ll[ . 
If our understanding is correct, the non-glassy regime 
is not necessarily fluidized, and the region of increas- 
ing granular temperature may overlap the region of non- 
vanishing configurational order associated with a "solid" 
state. This is an intriguing issue, and it is worth a sepa- 
rate study. 
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FIG. 4. Phase diagram of our system, and the damp- 
ing regimes, (a): explicitly analyzed regions, the value 
of the damping parameter b shown by grey scale vs the 
adimensional frequency and amplitude of shaking (white 
- highest damping). "Mountains" at low amplitudes are 
caused by fluctuations in the glassy state. The black 
stripe at the very top of the picture represents empty 
data and should not be taken into account. The diagram 
(b) is a summary of regimes based on the entire analysis. 
The narrow region between the glass state and dashed 
line is found to be a non-glassy "solid" . 

At the glass transition point, F « 1, particle ensem- 
bles get trapped in their momentary position. In the 
glassy regime the rate of dissipated energy is strongly 
sensitive to the configurational arrangement of the par- 
ticles i.e. to the history of the system. In other words, 
depending on initial conditions, for the same parameters 
of oscillation A and w, the system can settle in different 
configurational arrangements, and exhibit different dis- 
sipation rates. Transitions between different condensed 
phases based on multi-particle condensation have been 
suggested earlier by two of the authors ||] . 

To demonstrate this fact, in Fig. ^ two different sets of 
data points are shown by filled and open symbols; they 
refer to "fast" and "slow" cooling. By fast cooling we 



mean an instantaneous transition from any fluidized ini- 
tial state to the specified amplitude and frequency of vi- 
bration. The initial conditions, namely, random positions 
and velocities of particles, are kept the same for each data 
point. Thus, there is no memory of the preceding evolu- 
tion. Having the first 30 cycles of the driving oscillation 
discarded, the mean dissipation is obtained by averaging 
over 30 cycles of shaking. The large fluctuation of the 
dissipation rate b is not due to insufficient data averag- 
ing: averaging over 100 cycles, instead of 30, leads to 
the same results. In the case of slow cooling, the system 
is initialized only once for each frequency, and the final 
state for a certain amplitude A serves as the initial one 
for the next amplitude to be investigated. The ampli- 
tude of vibration in each simulation is diminished by a 
factor of 1.2. The behavior in this case can be sensitive 
to the entire history of the experiment. According to 
the mechanism discussed in Fig. in the region F ^ 1 
the fast cooling data points begin to fluctuate, whereas 
the slow cooling data are relatively smooth. Thus, for 
accelerations F ^ 1 one identifies a glass regime. 

The different damping regimes can be clearly appreci- 
ated in Fig. IJa, where a shaded contour plot of the ef- 
fective damping parameter is shown as a function of the 
adimensional amplitude and frequency of shaking [ p^ . 
These results suggest that, on the {A,uj) plane, the 
regime boundaries form a diagram similar to the Van der 
Waals s ystem , with the "triple point" roughly located 
at (L, y/gjL). Above this point, direct "sublimation" is 
achieved at Auj'^ « g. Convection rolls can develop at 
A < L above a critical velocity, as has been discussed. 
This diagram is shown schematically in Fig. where 
for the seek of clarity the boundaries between regimes 
are sharpened. The damping decrement along with con- 
vection dynamics provide indicative suggestions about 
"phases" which we discuss here. Any sound indentifi- 
cation of transitions and phases can only be based on 
studies of appropriate order parameters. In principle, 
it is conceivable that some "phases" may represent only 
preferred dynamic states. 

The fluidized regime is amendable to hydrodynamical 
arguments. Let us consider the "gas" state, where the 
dissipation is dominated by collisions with the walls of 
the container. The pressure transmitted to the vertical 
wall by the granular cloud of density p which moves with 
velocity v = Acu relative to the wall is ~ pA'^uj'^. Neglect- 
ing other contributions, the dissipated power during the 
collision is then roughly a fraction of the quantity pv'^L. 
Using Eq. (^) one finds 



(6) 



Formula (|^) is applicable at Au > a/C^. Cg is some 
unknown numerical prefactor. This expression is simply 
proportional to the collision frequency of independent gas 
particles moving at velocity Auj and having a mean free 
path of the order of L. In other words, collective modes 
are irrelevant in the "gas" phase as it should be. 
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For velocities Auj < \/Lg, the system is only partially 
fluidized. This means that the pressure transmitted by 
the container walls does not exceed the static stress that 
the condensed particles can sustain. The latter is roughly 
of the order of p — pgL in horizontal direction. As long 
as there is any fluidized material in the system, these 
stresses cannot be lower. Therefore, the dissipated power 
is pHv, where H is the height of the system. Dividing it 
by the total kinetic energy ^pHLv^ one gets 
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Auj 

Formula (0) is applicable at Aoj"^ > g, Auj > VRg, 
Auj ^ V^g, and C'l is the unknown numerical prefac- 
tor. Very roughly, in this regime some condensed par- 
ticles transmit their motion to other particles and the 
latter can move against gravity, b is the inverse time 
needed to decelerate them. The best fit of the curves of 
Fig. ^ provides the values Cg ~ 0.1 and Ci — 0.45 for 
our simulated system, giving a value for the minimum 
{Auj)a = 2.1y/Lg. This is obviously consistent with the 
upper curve shown in Fig. separating the "fluid" and 
"gas" regimes. 

To summarize, by means of MD simulations we have 
shown how the analysis of energy-loss rate displays dif- 
ferent damping regimes. In particular, one finds that 
fully convective states correspond to minimal damping 
decrement. Performing two types of measurements, re- 
ferred as to slow and fast cooling, one identifies a glass 
regime. This regime, in which configurational states af- 
fect the dynamical properties of the system, is separated 
from the fluidized regime by the value of the forcing pa- 
rameter, r = Auj'^/g ~ 1. For higher forcing one may 
find a non-glassy solid phase as long as the velocity of 
shaking is smaller than Auj < y/Rg- The Aw-scaling of 
the damping curves signals the beginning of the "fluid" 
regime. Here convective states can develop in a region of 
the plane (w. A) above the critical velocity Auj ~ VRg, 
critical acceleration Auj'^/g ^ 1 and below "evaporation" 
threshold A = \J Lgjuj. The damping decrement passes 
through the minimum for velocities Auj ^ ^/ Lg, and at 
higher velocities a gas-like state can be identified, where 
the effects of gravity are negligible. 
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